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Perturbative expansions of QCD observables in powers of a s are believed to be asymptotic and 
non-Borel summable due to the existence of singularities in the Borel plane (renormalons). This 
fact is connected with the factorization of scales (which is inherent to QCD and asymptotic free- 
dom) and jeopardizes the convergence of the perturbative expansion and the accurate determina- 
tion of power-suppressed corrections. This problem is more acute for physical systems composed 
by one or more heavy quarks. In lattice regulations, it reflects on the appearance of power-like 
divergences in the inverse of the lattice spacing for a series of quantities (A, gluelump masses, 
the singlet and hybrid potentials, ...) making that the continuum limit can not be reached for 
them. Nevertheless, all these problems are solved within the framework of effective field theo- 
ries with renormalon substraction. This allows us to obtain convergent perturbative series and to 
unambiguously define power corrections. In particular, one can connect with lattice results. Re- 
markably enough the dependence on the lattice spacing can be predicted by perturbation theory. 
This framework has been applied to the prediction of the gluelump masses and the singlet and 
octet (hybrid) potentials at short distances, as well as to their comparison with lattice simulations. 
Overall, very good agreement with data is obtained. 
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1. Introduction 

The perturbative series relating the pole and the MS mass: 

oo 

'"os =m M s + £r„(v)< +1 (v), (1.1) 

suffers from renormalon ambiguities, which makes this series asymptotic and non-Borel summable. 
The behavior of the perturbative expansion at large orders is dictated by the closest singularity to 
the origin of its Borel transform: 

B[m OS } (t(u)) = N « v ^_l u y +b (1 +ci(l — 2u)c 2 {\ - 2uf + •••) + (analytic term), (1.2) 
where {u = 



/oo ^ n 
dte- t/as B[m os }(t), B[i»os](0 = E r » — ■ (1-3) 
n n . 



o 



r, 

n=0 

The asymptotic behavior of the coefficients of the perturbative series then reads 

flr „^~ f PoY T{n+\+b) ( b bib -I) \ 

Quite remarkable, it is possible to obtain: 

a) the coefficients b, c\, cz,.., exactly through the use of the renormalization group [| | (actually 
they will depend on the coefficients of the beta function: /3o, j8i, ...). 

b) approximate determinations of the normalization constant, N m , ^ by defining new functions 
with improved analitical properties in the Borel plane [||], such that for those it is possible to 
perform an analytic expansion in the Borel parameter u. The determination of N m is quite solid and 
survives a series of checks, see 



• Good convergence of the perturbative series in u that determines N m . 

• Mild scale dependence of N m . 

• Consistency with the determination of Ny s , the normalization constant of the infrared renor- 
malon of static singlet potential V s . 2N m +Afy s — 0. 

• Agreement of the absolute value and scale dependence of the exact and asymptotic estimates 
of the coefficients of the perturbative series. See Fig. |]. 

2. Applications 

As we can see in Fig. [j], in heavy quark physics the asymptotic behavior sets in at quite low 
orders in perturbation theory, and the powers of ln[v/m] effectively exponentiate, becoming a linear 
power-like divergence in the factorization scale 1 : 



'At this stage the similitude with a 1/a power-like divergence one would find in lattice computations is evident. 
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r 2 v/m m r 3 v/m m 

Figure 1: Plots of the exact (rff) and asymptotic (r" s ) value of r„ (v) at different orders in perturbation theory 
as a function of v/m^g. The scale dependence of rf is known exactly. The constant term has been fixed 
using renormalon dominance. 



The associated lack of convergence of the perturbative series at low orders in perturbation theory 
becomes a problem in applications to heavy quark physics. The solution proposed in Ref. ^ was 
to shift the n\ factorial behavior from the perturbative series to the low energy matrix elements, 
where it properly belongs, since this behavior is associated to low energy dynamics. For the heavy 
quark mass, this implies to work with the RS mass: 

m RS (Vf)=m OS -8m RS =m os -'£N m V f ( j^Y a? +1 (v f ) £ c / ( "+ * + b ~ k) . (2.2) 

n =o \ l% J k=o J-U+o k ) 

This framework of renormalon subtraction can be applied to any effective theory with heavy quarks: 
HQET, NRQCD, pNRQCD, ... (see []|, [6|] for reviews). In this scheme some parameters become 
dependent on the scale, Vy , and scheme of renormalon subtraction. Here we will focus on a series 
of (quasi-) observables that can be studied in the static limit: 

(M B/D ) = m b/CjRS (v f ) + A RS (v / ) + 0(l/m b/e fls) , (2.3) 
E,(r) = 2m RS (v f ) + V,, RS (r; V/) + 0{r 2 ) V s = -C F a s /r +■■■, (2.4) 

E H {r) = 2m RS (v / )+K,Rs(r;v / )+A^ s (v / ) + ^(r 2 ) V a = \/{2N c )a s /r+ ■ ■ ■ , (2.5) 
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where H labels the hybrid/gluelump state at short distances, and the label RS for V s / , A, A# means 
that the leading infrared renormalon has been subtracted (added) from the perturbative series. The 
use of this scheme significantly improves the convergence of the perturbative series and the agree- 
ment with experiment and lattice simulations, when available, see [^|, |J [7|]. Examples include the 
determination of the IS bottomonium mass and the agreement with the singlet and hybrid poten- 
tials computed in the lattice at short distances. Another nice example is the accurate description of 
the dependence on the lattice spacing of the static singlet and hybrid potential, A and A# using per- 
turbation theory. This is possible because one may also consider to do the renormalon subtraction 
in a scheme different from the one presented here (as far as the subtracted quantity has the same 
non-analytic behavior in the Borel plane this would be legitimate). Here we would like to connect 
with the lattice scheme. In practice, this means to perform the replacements: {RS — > L, Vf — > 1/a}. 
The relation between both schemes is renormalon free and governed by a convergent series as we 
shall next see. 



3. Perturbative running of 1 /a and scheme dependence 

Due to the existence of renormalon ambiguities (1/a divergences in the lattice language), it is 
not possible to get the continuum limit of the static singlet and octet potential, gluelump masses, 
Ah and A. Nevertheless, their dependence on I /a is free of renormalons and can be predicted by 
perturbation theory: 

2(A L (l/a) - A L (l/a')) = V SiL (r,l/a) - V s r(r;l/a') = C F ( - - - f ) Vl a s + ■ ■ ■ . (3.1) 



The coefficients vi, v%, V3 are known in pure gauge theory with Wilson action [Jsj, ^J. A similar 
renormalon-free equality can be constructed for the octet potential and gluelump masses. 

A H {\/a)-A H {\/a') = [V 0jL (r,l/a) -V s>h (r;\/a)] - \V 0>h {r;\/a!) -V s>L (r,l/a')] 

C A (\ 1 



, , 7 W + (3.2) 
2 \a a! ) 

where v \ , v' 2 are exactly known and V3 in the large N c limit. 

We show the plot corresponding to Eq. ( 3T ) in Fig. ||, where the static potential data has been 



normalized to agree with A at one specific value of the lattice spacing 

Apot(«) = ^ L (^o;a)+A. (3.3) 



We can see how nicely the perturbative prediction in the lattice scheme (continuous green line in the 
first plot in Fig. 2) agrees with the lattice data [[R], §,11]. Significantly, only for one collaboration 



the slope of the lattice data points slightly differs from the prediction of perturbation theory [gj. We 
can also see the convergence of the perturbative series relating the lattice and the RS scheme, since 
it is also renormalon free. The final value for Ars, Ars(V/ = 2.5r^ l ) = 1.17rQ , agrees within 
errors with the result obtained directly from experiment using a combined analysis of the T(15) 
and the B meson mass [|, A RS (v/ = 2.5r l ) = 0.92r l . 

A similar analysis can be performed for the octet potential and gluelump masses []7|]. In this 
case there is less statistics but the slope predicted by perturbation theory nicely agrees with the 
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Figure 2: The first figure shows the binding energy Ap 0t (full diamonds [0), in comparison with A 

(open diamonds ||, 11 1). NLO, NNLO and NNNLO refer to transformations of Ap Gt into the RS scheme to 
different orders in perturbation theory. The solid line corresponds to the NNNLO evaluation in the lattice 
scheme with the central value A (V/ = 9.76r ( ^ 1 ) = l.VQrjj" 1 . The second figures shows A pot translated into 
the RS scheme at NLO (squares), NNLO (pentagons) and NNNLO (circles). The solid line corresponds to 
the principal value running in the RS scheme. The error band corresponds to the prediction A (9.76 r^ ) = 
(1.70 ±0.04) rg 1 , and includes the uncertainty due to A m = (0.602 ± 0.048) r Q l . For details see [^|. 
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one obtained from the gluelump mass A| measured at different lattice spacings [12]. Moreover, 
the perturbative series governing the change to the RS scheme converges well. Finally, it is quite 
comforting that the same value, within errors, is obtained for Ag from either the hybrid potential, 
A| s (v/ = 2.5r 1 ) = 2.25r 1 , or from the direct lattice determination, A| s (v/ = 2.5r 1 ) = 2.3 lr Q 1 . 
All these findings can be found summarized in Fig. 14 in [|7|]. 

In brief, we can relate processes computed with different scales/schemes using well-behaved 
(renormalon free) perurbative series: 

A L (l/a) Ai(l/«') 1 

§. (Vf = 1/a) § (v'f = 1 /a) > The circle can be closed using perturbation theory 

Ars(vj) Ars(v}) J 

and a similar circle applies for V s L (r; 1 /a), (r; 1 /a) and A^(l/a). 



4. Conclusions 

We have accumulated a lot of evidence in favour of the renormalon dominance in heavy quark 
physics. A proper handle of these effects appears to be crucial to accurately describe either lattice 
or experiment. We point out that the dependence on the lattice spacing can be obtained from 
perturbation theory with good accuracy and well controlled errors. Therefore, lattice simulations 
in heavy quark physics can be peformed with quite coarse lattices and yet obtain accurate results. 
This may have important consequences to diminish errors in lattice simulations. 
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